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Abstract 

It is well-known that Noether currents in the classical four- dimensional M = 
1 supersymmetric Yang-Mills theory (4D M = 1 SYM), i.e., the U(1)a cur- 
rent, the supersymmetry (SUSY) current and the energy-momentum tensor, 
form a multiplet under SUSY, called the Ferrara-Zumino supermultiplet. In- 
spired by this structure, we define the energy-momentum tensor in the lattice 
formulation of 4D M = 1 SYM by a renormalized super transformation of a 
lattice SUSY current. By using a renormalized SUSY Ward-Takahashi rela- 
tion, the energy-momentum tensor so constructed is shown to be conserved 
in the quantum continuum limit. Our construction of the energy-momentum 
tensor is very explicit and usable in non-perturbative numerical simulations. 

Keywords: Lattice gauge theory, Supersymmetry, Energy-momentum 
tensor 



1. Introduction 

Ideally, a non-perturbative formulation of a field theory with some sym- 
metries should provide, not only the definition of correlation functions, but 
also the definition of renormalized Noether currents that generate correctly- 
normalized symmetry transformations on renormalized fields. This is ex- 
pressed by renormalized Ward-Takahashi (WT) relations and, when these 
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relations hold, one may say that the symmetries are really realized in quan- 
tum field theory. 

Quite often, however, the regularization procedure breaks the preferred 
symmetries and, for this reason, it is generally very difficult to conclude even 
the existence of such renormalized Noether currents (especially when the 
symmetry transformations are non-linear). Even if such Noether currents 
are assumed to exist, the explicit construction can be very cumbersome. In 
the lattice formulation of supersymmetric theories, one encounters such a 
situation, because the infinitesimal translation (that is a part of the SUSY 
algebra) is broken by the spacetime lattice, implying also the breaking of 
SUSY. Here, almost all fundamental symmetries that define the theory are 
broken by the regularization. 

Having the above general remark in mind, in the present paper, we study 
the construction of Noether currents in the lattice formulation of the four- 
dimensional M = 1 supersymmetric Yang-Mills theory (4D J\f = 1 SYM) [l[ 
(see also Ref. j^j for an earlier consideration and Ref. p| for a very readable 
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review). The classical Euclidean action of 4D N = 1 SYM is given by0 

^1™! = j d A x ji tr [F^{x)F^{x)\ + tr $(x)pij>{x)] | , (1.4) 

where the adjoint fermion (gluino) ip(x) on the Euclidean space is subject of 
the constraint 

^(x)=^ T (x)(-C" 1 ), (1.5) 

to express the degrees of freedom of a Majorana fermion in the 4D Minkowski 
space. The global SUSY 5 C in the classical 4D A/" = 1 SYM is 

1 - - 1- 

8^{x) = --cr^F^(x), S^(x) = -fr^F^lx), (1.6) 

where the Grassmann-odd constant spinor £ obeys £ = £ T (— C -1 ). Without 
containing the scalar field, 4D M = 1 SYM is the simplest supersymmetric 



1 We basically follow the notation of Ref. [4|: The sum over repeated indices is under- 
stood. Vector indices fi, v, . . . , run over 0, 1, 2, 3. e^pa- denotes the totally anti-symmetric 
tensor and £0123 = — 1. All gamma matrices are hermitian and obey {7^,7^} = 2<$ M „. We 
define 75 = —70717273 and a^ Ll/ = [y fl ,j l ,]/2. The charge conjugation matrix C satisfies, 
C _1 7 M C = -7J, C^a^C = -o~£,, C~ l ^C = 7 J and C T = -G. The generator of 
the gauge group SU(N C ), T a , is normalized as tr(T Q T b ) = (l/2)S ab . F^{x) is the field 
strength F llu {x) = d^A v {x) — d v Ap{x) + ig[A /1 (x), A„(x)], where g is the bare gauge cou- 
pling constant, and p = 7 P -D M , where = ig[A^{x),] is the covariant derivative 
in the adjoint representation. For the lattice theory, x, y, z, . . . denote lattice points and 
a is the lattice spacing; fi is the unit vector in the /x-direction. The link variable U^(x) 
and the gauge potential A^x) are related as 

U ll {x) = e iaBA ^ x \ (1.1) 

The forward and backward difference operators respectively are defined by 

d^f(x) ee I [f(x + afi) - f{x)] , d;f(x) ee 1 [f(x) - f{x - ay)] , (1.2) 

and the symmetric difference operator is defined by 

d s ^EE l -id, + d;). (1.3) 
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theory in four dimensions and we should understand the issue raised at the 
beginning of this paper first in this example, before tackling more complicated 
supersymmetric theories. 

The most interesting result we will obtain below is a very explicit form 
of an energy-momentum tensor on the lattice T pv {x), given by Eqs. (13 . 3[) 
and f)3.4p . which is conserved in the quantum continuum limit. Our way 
of construction of this energy-momentum tensor was suggested by the fact 
that in the classical 4D M = 1 SYM, the Noether currents associated with 
the U(1)a symmetry, SUSY and translational invariance, form a multiplet 
under SUSY, the Ferrara-Zumino (FZ) supermultiplet For the classical 
theory ( II. 4p . we define the U{1)a current %^{x), the SUSY current S p (x) 
and the energy- momentum tensor T pu (x) by 

j 5p (x) = tr [^(2)7^75^(2;)] , 
Sf,(x) = -o-p^tr [ip(x)F pa (x)) , 

f pu (x) = 2tr [F pp (x)F up (x)} - -(^tr [F pa (x)F pa (x)} 

ip(x) \^f p D U + ^ U D p Jip{x) --5 pu tr ip(x) Tp 



1 

+ - tr 
4 



(1.7) 



where the definitions of j5^(x) and S p (x) are standard while, as well-known, 
the definition of the energy-momentum tensor is large extent arbitrary. In 
Eq. (II. 7p . we defined the energy- momentum tensor by first introducing a 
background gravitational field into the continuum action (jl.4p and then tak- 
ing a flat-space limit after differentiating the action with respect to the grav- 
itational field. Then one finds, under the super transformation (11. 6p . 

f>t35»{x) = €7sSp(x), (1.8) 



2 See Chapter 20 of Ref. [6| for a very readable exposition. An interesting application 
of the notion of the FZ multiplet in phenomenology was recently considered in Ref. 0] ■ 
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and 



+ Ihlvidvjs^x) 

3 3 - 

- -xlhlui tr l^(x)^^^ u pi)(x)] 

- |0V P 6tr [^(x)^^a up Ipi)(x)\ . 



1.9) 



Similarly, 5{T pu {x) becomes a linear combination of d p S a (x), up to terms be- 
ing proportional to the equation of motion of the (massless) gluino, pip(x). 
By integrating over the spatial coordinates of the \x = component of the 
above relations, one sees that they are consistent with the defining algebra 
between the SUSY charge, the U(1)a charge (this is the i?-charge) and the 
momentum, up to the equation of motion pip(x) = 0. Eqs. (jl.8p and (ll.9p 
thus may be regarded as a basic characterization of fundamental Noether 
currents in Eq. (|1.7p . By going one step further, one may regard Eqs. (I1.8P 
and (II. 9p as defining relations of the Noether currents. This is the approach 
we adopt in the present paper. That is, we define a lattice energy-momentum 
tensor T^vix) by a (renormalized modified) super transformation of a (renor- 
malized) lattice SUSY current S^ix) that is conserved in the continuum limit 
under an appropriate gluino mass tuning [lj. Then it can be shown that the 
lattice energy-momentum tensor so constructed is conserved in the contin- 
uum limit. This conservation of the energy-momentum tensor is not quite 
trivial, because SUSY is not a manifest symmetry with the lattice regu- 
larization. For the conservation of 7^^(x), a "finiteness" of a lattice super 
transformation induced by S^ix) turns to be crucial. This finiteness might 
be regarded as a necessary condition for the existence of a renormalized 4D 
M = 1 SYM. 

As the structure of our energy- momentum tensor is very explicit as Eqs. (13. 3p 



and (13.41) s 



SYM 



23,1221 



row, it is usable in actual Monte Carlo simulations of 4D M = 1 



a]9|. 

,13, 




J? 



nL H HHEE E3, EE EU S 0, H S H [22, 

32[. Physical quantities such as the viscosity (see, for exam- 
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pie, Ref. 33J and references therein) might be measured from the two-point 
function of our lattice energy-momentum tensor. 

Our argument below is quite general and rather independent of the de- 
tails of a lattice formulation. We assume the locality of the lattice action, 
that consists of the sum of a gauge boson (gluon) action Sgiuon, the kinetic 
term Sgaho and the mass term S masB of the gauge fermion (gluino) Jf| 

S — S'gluon 'S'gluino 'S'mass; (^■"^■^) 

where 

Sgiuino = a 4 ^tr $(x)Di(>(xj\ , (1.11) 

x 

where D is a lattice Dirac operator ((C _1 D) T = — C _1 .D), and 

<? mass = a 4 ^Mtr ${x)if>(x)] . (1.12) 

X 

We further assume that the lattice action fll.lOp is invariant under the hy- 
percubic group H(A), including the parity transformation V defined by 

U (x , x) —> U (x , —x), U k (x , x) —> Ul(x , —x - ak), 

t/j(x , x) -)> ij if>(x , -x), i>(x , %) -> -iip{x , -0)70, (1.13) 

where and in what follows, x = (xi, 22, £3) and k, I = 1 or 2 or 3. 

In the following discussion, we make frequent use of a lattice transcription 
of the field strength, 

[F, U } L (x) = 2tr [P, u (x)T a ] T\ (1.14) 
that is defined from the clover plaquette P^x), 

1 4 1 r i 

P ^) E lEFi- U^-UlM , (1.15) 
4 L — ' Liarq L p J 



3 Although we use the language of a four-dimensional lattice theory in the present 
paper, we believe that our argument can be transcribed for a domain-wall-type [34], |35| 
five-dimensional setting 3^, 37 , 38 1 without much difficulty. 
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where 

U lflu (x) = U^(x)U v (x + ap)Ul(x + au)Ul(x), 

U2nv(x) = U v (x)WJx — a/t + aO)Ul(x — ajl)U^(x — a/i), 

U3fi V (x) = W{x — ap)Ul(x — aft — au)U^{x — aft — au)U v (x — az>), 

U^ u {x) = Uj,(x — au)U^(x — av)U v {x + a/t — av)W{x). (1-16) 

Note that [F AlI/ ] I '(x) is traceless by construction (I1.14p . Under the parity, 
[F/iv} L (%) transforms in an identical manner as the continuum field strength. 
That is, we have 

[F ok } L (x , x ) A - [F Qk } L (x Q , -x), [F m } l (x , x) A + [F H } L (x , -x). 

(1.17) 



2. Lattice SUSY current and a renormalized SUSY WT relation 



In this section, we recall how the proposal of Ref. [1| that SUSY in the 
lattice formulation of 4D M = 1 SYM is restored with an appropriate tuning 
of the ghiino mass M is understood in terms of the SUSY WT relation |l|, 



39, 20 



See also Refs. 17, 19 



This discussion provides the definition 
of a lattice SUSY current that is conserved in the continuum limit and our 
basic SUSY WT relation that will be fully utilized in the construction of the 
energy-momentum tensor in Sec. [3j 

First, as a lattice counterpart of the SUSY transformation in the contin- 
uum (11. 6p . we define |39 



S^U^x) 
8^{x) 



1 - 



(x)Un(x) + U^(x)ijj(x + a/t)] 



= a 



[X 



5&{x) 



[X 



(2.1) 



To derive a WT relation, we need also the localized transformation and, 
corresponding to the above, we define 

S^U^x) = iag^ [l(x)-i^(x)U^(x) + l(x + ajx)^JJ^(x)il)(x + ajx)] , 
s ( u t( x ) = - ia d\ \^{ x )ln U l{x)il){x) + & + aflj^x + afi)Ul(x)] , 



8^{x) = --tr^{x) [F^] 1 (x), 5^(x) 



\l{x)a lxv [F, v ] L {x), (2.2) 
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where the Grassmann-odd spinor parameter £(x) obeys £(x) = £ T (x)(— C^ 1 ). 

Now the lattice action (11 . 10p is not invariant under 5% (12. ip . because, first 
of all, the gluino mass term (jl . 12j) explicitly breaks SUSY and, various 0(a) 
lattice artifacts in the lattice action also break SUSY. Since the localized 
SUSY transformation 5^ (12. 2p reduces to the global one (12. ip for — > £, 
the variation of the lattice action (ll.lOp under Eq. ( 12. 2 p can be written as 

6^S = a 4 ^£» [-tfS^x) + M x (x) + X s (x)] , (2.3) 

X 

where combinations S^(x) and x( x ) are defined by 

S^x) = -a pa%l tr fax) [F p(7 } L (x)} , (2.4) 

and 

X {x) = a^txfax) [i^] L (x)} , (2.5) 

respectively. S^x) in Eq. ( 12 .4p is nothing but a lattice transcription of the 
SUSY current in the continuum, S p (x) in Eq. (II. 7p . and x( x ) i n Eq. (12. 5p is 
an explicit SUSY breaking caused by S mass . Thus, in Eq. ( 12. 3p . X s (x) repre- 
sents an 0(a) SUSY breaking effect attributed to the lattice artifacts. Thus, 
considering the infinitesimal variation (I2.2p in the integration variables of the 
functional integral, we have an exact identity that holds for any operator 

(d*S»(x)0) = (\M X (x) + Xs(x)\ O) - (^J*8(Pj . (2.6) 

Here is a remark on the property of Xg(x) under the parity: From 
Eqs. and ([23]), we see that 

o / — *\ *P (+hoS (x , -x) for fj. = 0, 
S„(x ,x)^< (2.7) 
{-HoSk(x , -x) for /i = k, 

and 

X(x , x) ->■ hox(xo, ~x). (2.8) 



4 Throughout this article, we always assume that an operator represented by the sym- 
bol O is gauge invariant. 
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These shows that S^(x) possesses the same transformation property as x( x ) 



under the parity. On the other hand, by defining ^(xq,x) — > i r y ^(x ,—x) 



■p 

and £(xo,x) — >■ — i£,(xo, — x)7o, we see that ^ in Eq. ( 12. 2 p and the parity 
commute to each other. Then by the parity invariance of the action S and 
Eqs. (EZD and in Eq. Q obeys 

X 5 (x , z) ->■ i7 X 5 (x , -f). (2.9) 

Next, we consider the renormalization of the composite operator Xg(x) 
(just like for the chiral symmetry 40, 41]). It can be argued that [H,0,0|, by 



using the assumed symmetries of S and resultant properties of Xs(x) (such 
as Eq. O), by writing^ 

X s (x) = (1 - Z s )d s ^(x) - Z T d%{x) - -Z xX (x) 

(X 

- Z 3F ti [ip(x)i)(x)ip(x)] 

- Z EO M^tr{[F^] L (x)( J D + M)iP(x)} 

+ a£(x), (2.10) 

and choosing renormalization constants Z s , Z T , Z x , Z 3F and Z EOM appro- 
priately, the dimension 11/2 operator S(x) can be made at most logarithmi- 
cally divergent for a — > 0. More precisely, the operator £(x) is given by a 
linear combination of renormalized operators with logarithmically divergent 
coefficients. Since the operator O in Eq. (12.61) is gauge invariant, the mixing 
of Xs(x) with gauge non- invariant operators [39( does not occur in Eq. (I2.10p . 
A new lattice operator T M (x) in Eq. (I2.10p is defined by 

T M (x) = 2 7v tr [iP(x) [F, U ] L (x)} , (2.11) 



5 In Ref. [J], this operator renormalization was described by using an operator basis 
in the continuum theory. With the present lattice regularization, however, it is more 
appropriate to use an operator basis in the lattice theory as Eq. (|2.10[) , because the renor- 
malization is carried out in terms of counterterms with lattice structure [42l l43l l44j . Also, 
we adopt a different choice of basis for the term proportional to the equation of motion 
from Ref. as this choice is practically more convenient. Finally, we make use of a 
simpler lattice SUSY current than Ref. [39(. Then the symmetric difference instead of 
the backward difference d* appears for the covariance under the parity. 
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and, from this definition, 



m , -« v K J +^7o^o(a;o, -x) for /i = 0, 
[-z7 T fe (x , -ac) for // = fc. 

From Eqs. ( 12. 6 p and (I2.10p . we have 

(d s ^[Z s S,(x) + Z T T,(x)}0) 
1 
o 

1 d 



(2.12) 



M - -Z x ) (x(x)O) - Z 3F <tr [^(x)^(x)^(a;)] 0) 

<^cA - Z EOM (<V ti{[F^] L {x){D + M)i/;{x)}0) 



a 4 <9£(z 
+ (o5(ar)0) . 

In this expression, we first note 



a 



[MP 



tr{[F^] L (x)(D + M)^(x)}0) 







a A d£(x) 



(2.13) 



(2.14) 



where the localized transformation 8 F £ is a super-like transformation that 
acts non-trivially only on the gluino: 



Sf^U^x) = 0, 5 F £ip(x) = <%0(z), S F ^ip(x) = 5gil>(x), 



(2.15) 



because the left-hand side of Eq. f |2.14j) is proportional to the equation of mo- 
tion of the gluino (i.e., this is the Schwinger-Dyson (SD) equation). There- 
fore, we have 

(dS[Z s S,(x) + Z T T,(x)}0) 

M - -Z x ) ( X (x)0) - Z 3F <tr ty{x)${x)il>(x)] O) 







-Af + aS(x) 



O 



a 4 d£(x) 

where A^ is a modified localized super transformation, 

A^ = 5^ + 2eom£f£- 



(2.16) 



(2.17) 



Let us consider a special case of Eq. (I2.16P that O is a collection of 
renormalized local operators and the point x stays away from the support 
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of O by a finite physical distance. In what follows, we express this situation 
by 

x *vv» supp(O). (2.18) 

Under this situation, the £ (^-derivative in the last line of Eq. ( I2.16P identi- 
cally vanishes^) Also the last term of Eq. (I2.16P vanishes in the continuum 
limit, because for x supp(O) the dimension 11/2 operator £(x) does not 
produce any 0(1 /a) divergence that can compensate the overall factor of a. 
We thus have 



(d s ^[Z s S,(x) + Z T T,(x)]0) 

^ (m - -Z x ) (x(x)O) - Z 3F <tr ty(x)$(x)^(x)] O) 



a 

for x supp(C). (2.19) 

This can be regarded as a would-be conservation law of a lattice SUSY cur- 
rent. If this relation with Zsf ^ held, then the last term, that is cubic in 
the gluino field, would give rise to an "exotic" SUSY anomaly. This is not 
what we expect and we thus assume its absence: 

Z 3F = 0. (2.20) 

As shown in Ref. [If by utilizing the generalized BRS transformation, this is 
actually the case at least to all orders in the perturbation theory. Accepting 
Eq. (127201) . we then have 

[Z s S,(x) + Z T T,(x)] O) ^ (m - -Z^ (x(x)O) , 

for x «v«* supp(O). (2.21) 

The combination M — (l/a)Z x in the right-hand side corresponds to an 
additive renormalization of the gluino mass. So we tune the bare mass pa- 
rameter M [1] so that 

M--Z x = 0. (2.22) 

a 



6 Our argument below works when the operator O is ultra-local, i.e., its support is 
a strictly-finite region on the lattice. With the use of the overlap lattice Dirac opera- 



tor D [45|, |46|, however, one might be interested in an operator O that is not ultra- local 
but exponentially local [47j . Although we do not give any analysis for such a case, on 
physical grounds, we believe that our conclusions will not change. 
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In actual numerical simulations, this tuning should be carried out for each 
value of the lattice spacing a. In this way, we obtain the conservation law in 
the continuum limit: 

(9J \Z s S,(x) + Z T T^(x)] O) ^ 0, for x <— supp(O), (2.23) 

that is a minimal requirement for the restoration of SUSY in the continuum 
limit. 

We consider the renormalization of the conserved current ZsS^(x) + 
Z T T^(x). Since this is a dimension 7/2 gauge invariant vectorial operator, 
it can mix only with the operators S^x) and T^(x). However, if the mixing 
were generic, it would be inconsistent with the conservation law ( I2.23p . That 
is, the renormalization must be multiplicativeJll Thus we set 

<S M (x) = Z [ZsSpix) + ZtT^x)} , (2.24) 

and call this a renormalized lattice SUSY current. With an appropriate 
choice of the constant Z, S^x) should have finite correlation functions with 
any renormalized operator O, as far as x <^ supp(O). The renormaliza- 
tion constants Z$ and Zt in Eqs. (I2.10p and fj2. 24j) emerge from the power- 
divergence subtraction and it can be argued that such constants are indepen- 
dent of the renormalization scale fi 2 and finite as a — > [41] . On the other 
hand, by the dimensional reason the constant Z in Eq. (I2.24p is at most log- 
arithmically divergent. In terms of the renormalized SUSY current (I2.24p . 
the SUSY WT relation fl2~T6|) under the assumptions fT2~20p and fl2T22|) reads 

(d*S»(x)0) = (z 

Finally, we argue that for any renormalized local operator O, the combi- 
nation 

O, (2.26) 

when x G supp(O), produces a renormalized local operator in the support 
of O. In Eq. f !2.26p . something non-trivial in the continuum limit can occur 
only when x G supp(O), because otherwise the ^(x)-derivative vanishes and 



7 It is very interesting to confirm this renormalizability by the lattice perturbation 
theory along the line of Ref . [33] • 



1 d 

"a^Wx 



-A f + a£(x) 



O 



(2.25) 



Z 



d 



a 4 d£(x 



-Af + a£(x) 
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the dimension 11/2 operator S(x) does not produce any 0(1 /a) divergence 
that can compensate the overall factor of a. To see what happens when 
x G supp(O), we sum Eq. (12.251) over the point x within a finite region T>o 
that contains the operator O. Then, since EigPo^^W lY1 ^he left-hand 
side becomes a collection of (mutually non-overlapping) renormalized opera- 
tors that have no overlap with O, the left-hand side remains finite as a — > 0. 
This shows that the sum ^2 X&V °^ ^ ne right-hand side of Eq. (12.251) is finite 
too and we see that the combination (12.261) does not produce any ultravi- 
olet divergence even for x G supp(O). We note that such a finiteness will 
be indispensable for the existence of a renormalized SYM, because for that 
there should exist a renormalized SUSY current that generates a finite super 
transformation on renormalized fields. Since the above analysis shows that 
Sfj,(x) is the unique conserved SUSY current in the present framework, S^ix) 
should generate such a finite super transformation. This argument on the 
basis of the existence of a renormalized SYM also supports the finiteness of 
the combination (I2.26p . 



3. A lattice energy-momentum tensor and its conservation law 

We now define an energy- momentum tensor T^ v (x) on the lattice through 
the relation quite analogous to Eq. (II. 9p . That is, 

ZA^S^x) = 2j v £{T ltt ,(x) + cS^tr [$(x)(D + M)^(x)} 

+ (terms anti-symmetric in \i and z/)j 
+ (terms proportional to 757^ £, 7s£, 0VpO> (3- 1 ) 

where c is a constant and Ag is a "global version" of the localized super 
transformation (I2.17p . i.e., 

A ? ee 5$ + ZeomSf^, (3.2) 

and 5f£ is obtained by setting £(x) — > £ in Eq. (I2.15p . One can invert 
Eq. (13. ip with respect to 7^, (2): We introduce (here and in what follows, 
indices a, (3, . . . , which run over 1, 2, 3, 4, denote the spinor indices), 



1 d 

0^(x) ee -( 7v ) pa — [ZA^(x)] a 



Z 2 Z s -( y -f u )f3 a -^- I fa + Z EOM 5 F t) 



S^x) + (x) 
Zs 



(3.3) 



13 



where in the second equality we have substituted Eqs. (13. 2p and (I2.24p . then 



%„(x) = - [6V(x) + G Vfl (x)] - c^tr [$(x)(D + M)^j(x)} . (3.4) 



From relations obtained so far, one can confirm that O^ix) in Eq. (13.31) 
behaves as a second-rank tensor under the parity!^ 



■p 



+6>oo (Xq, -x) 

-O k0 (x , -x) 
-O ok (x , -x) 
,+O k i(x , -x) 



for 


n = 


u = Q, 






for 


» = 


k and 


V 


= 0, 


for 


V = 


and 


V 


= k, 


for 


V = 


k and 


V 


= I. 



(3.5) 



The defining relation (13. ip is quite analogous to the classical relation ( II. 9p . 
Yet they differ in that the super transformation 8^ is replaced by the renor- 
malized modified super transformation ZA^ and the Dirac operator is shifted 
by the gluino mass M; also the coefficient of the term that is proportional 
to the fermion action is changed from 3/4 to c. 

Quite interestingly, we can show that the energy-momentum tensor de- 
fined by Eq. ( 13. 4p is conserved in the quantum continuum limit. That is 



(d% u (x)0) ^ 0, for x <~» supp(O), 



(3.6) 



for any renormalized local operator O. Note that this conservation law holds 
for any value of c in Eq. ( 13. 4p because the term tr[ip(x)(D + M)ip(x)} is 
proportional to the equation of motion and it has no correlation with O 
when x is not in the support of O. When x ^ supp(O), one can replace it by 



<tr [ip(x)(D + M)ip(x)] ) = -2(N 2 C - l)a 



(3.7) 



Since this is independent of x, the last term in Eq. (13 .4p does not affect 
the conservation law. The most general structure of the energy-momentum 
tensor (that is symmetric in its indices) in the lattice formulation of the 
Yang-Mills theory coupled to fermions has been given in Ref. 48] ; see also 
Refs. (49I . 50, 51 1 for related analyses. According to Ref. 48], a symmetric 



8 To see this, it is convenient to note that and the parity commute to each other, if 

V - v - 

one assigns £ — > 170^ and £ — > — i?7o- 
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energy-momentum tensor that satisfies the conservation law in the continuum 
limit ( 13. 6p . if it exists, is essentially unique. Only ambiguities are the overall 
normalization (Z 2 Zs in Eq. (13 .3p in our construction) and a proportionality 
constant to the fermion Lagrangian (i.e., the constant c in Eq. (|3.4|1 ). The 
bottom line is that our definition (13 .4p . that is suggested from the structure 
of the FZ supermultiplet, gives rise to a very explicit form of this unique 
conserved symmetric energy-momentum tensor. 

To show the conservation law (13. 6p . we first show the conservation law 
of 0^v(x) in Eq. (13. 3p that is not necessary symmetric under /if>i/. Then we 
show the anti-symmetric part of <9 MJ ,(x) also is conserved in the continuum 
limit, implying Eq. (13. 6p . 

To show the conservation law of 6> MJy (x), we set O —> d^S u (y)0 in the 



The second equality holds because the first expression is anti-symmetric un- 
der the exchange, x •<->■ y, //■<->■ v and a •<->■ j5. Using Eq. (I2.25P once again in 
the last line, we have the identity 



WT relation f[2^5j) . This yields 



([#)UM^) 




(3.8) 





a 4 d£(x) 



1 d 



A 5 + a£(x) Z 



« 4 a&y) 



1 d 



A f + a£(y) O 
- P 



(3.9) 
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We then sum this relation over y within a finite region T> x containing the 
operator d^S^x). Noting the identity between the local and global transfor- 
mations, 



yev. 



dtp 



(3.10) 



we have 







Z 



1 d 



a 4 d£(x 
1 d 



A e + a£{x) 



a 4 ]T [d s v SM],0 



Z 



a A d^{x 



-A^ + a£{x 



i 



a yev x 



1 d 



\ + aS{y) 



(3.11) 



Since £ = £ T (— C x ) implies 
d 



(7 







dtp 



dip 



x 



(3.12) 



Eqs. (E3D and f l3TTj) yield, 

<«9^(x)0) 



-(c-^u/z 



1 ^ A 

"1 a?7IX 5 + a£(z) 



--(C-V)aM-2 



a 4 d£(a; 
1 d 



a yev x 



a A d£(x 
yeV x 



1 g 



A 5 + a£(y) 



O 



(3.13) 



Now, let us suppose that the point x stays away from the support of the 
renormalized operator O by a finite physical distance, i.e., x supp(C). 
Let us further assume that the region T> x is chosen so that it does not overlap 
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with the support of O. In this situation, since ^2 y£T>x d^S u (y) does not have 
any support at x, Eq. (13. 131) reduces to 



\{C- l lv U(z[a£{x)] a a* £ [d s v S v (y)]^0 

d 



1/3 

y&Dx 




+ a£(x) 



Z[a£{y)] O 



a 4 d£{x) 

(3.14) 

The first term in the right-hand side is a correlation function of renormalized 
operators without any mutual overlap with the overall factor of a. Thus, 
this term vanishes in the a — > limit. In the second term in the right- 
hand side, we use the property argued around Eq. (I2.26p . Then it is also a 
correlation function of renormalized operators without any mutual overlap 
with the overall factor of a. Thus this also vanishes in the a — > limit. In 
this way, we conclude 

{dS0^(x)O) ^ 0, for x supp(C), (3.15) 

for any renormalized operator O. 

Next, we consider the anti-symmetric part of O^ix), 

A^(x) = ~ [0^(x) - e ni {x)} , (3.16) 

and its renormalization. It must be possible to expand this composite opera- 
tor by gauge invariant operators with dimensions less than or equal to 4 that 
behave in the same way under the hypercubic group as the right-hand side 
of Eq. (I3.16P (for the parity, see Eq. (|3.5|) ). By taking also the constraint (II. 5p 
into account, the most general possibility turns to be^| 

*A[iv (■£) 

= A 1 e tivpa d^ tr [ip(x)w 5 ip(x)] + A 2 tr [i/j(x)a^(D + M)^{x)] + aG^(x). 

(3.17) 



9 For example, a seemingly-obvious candidate ti[ip(x)a fi ^'ip(x)} identically vanishes be- 
cause of Eq. (|1.5|) : another candidate tr['tp(x)^ l j l D^tp(x)] — (fj, O v), where D^f(x) = 
(l/2a)[U fl (x)f(x + afyU^x) — U fJl (x — a/t) 1 ' f(x — afi)U^(x — a/t)], reduces to Eq. (|3.17jl . 
by the fact that D = ft + O(a) and a^ft = -e^p^Dp'j^ + J^D U - 7„-D M . 
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That is, by choosing the constants A% and A 2 appropriately, the dimension 5 
operator Q^ v (x) can be made at most logarithmically divergent. Now, quite 
fortunately, Eq. (I3.17P shows that the anti-symmetric part is conserved by 
itself: 

(d^A^(x)O) ^ 0, for x supp(O). (3.18) 

This is trivially true for the first term of Eq. (I3.17p . For the second term, this 
follows from the SD equation (the equation of motion of ip(x)). Finally, the 
last term of Eq. ( 13 . 1 7[) vanishes in the continuum limit. Combined Eqs. ( 13 . 1 8[) 
and f )3.15p . we conclude that the symmetric part of 0^ v (x), (l/2)[0^ u (x) + 
UfJi (x)}, also is conserved. This proves our assertion Eq. (13. 6p . because the 
symmetric part of O^x) is the first term of Eq. (\3A\i and the term with the 
coefficient c does not affect the conservation law as we already explained. 

We note that, with our definition of the energy-momentum tensor (13 .4p . 
the trace anomaly 52l. |53|. 154. |55| and the gam ma-trace anomaly (supercon- 
formal anomaly) [Ha" 57, 58, 59, 60, 61, 62, 63, 64] are related b}o 



(T m (x)0) 



d_ 

8 \dla 



[ZA^^S^ 



x)\ a O 



+ 8c(iV c 2 - l)a~ 4 (O) , (3.19) 



where we have used Eq. (13. 7p assuming x supp(C). 

Now we discuss how our energy-momentum tensor (13 .4p can be used in 
actual non-perturbative Monte Carlo simulations. As Eqs. (13. 3 p and (13.41) 
show, the definition of the energy-momentum tensor contains four combi- 
nations of renormalization constants, Z 2 Z$, Zeom, Zt/ ' Z$ and c; we have 
to know these numbers to construct the energy- momentum tensor. Among 
these, the ratio Z T j Zg has been measured [66j, |l6j, [l8|, |22j, |25J by determining 
the coefficients in the relation (O is a fermionic spinorial operator) 



d s ^s,(x)0} + ^ (y;,(.r)(0) 



1 



M--Z x \ (x(x)O) 



(3.20) 



This is Eq. (I2.2ip up to 0(a) corrections. Thus we already know that the 
number Zt/Z$ can be determined numerically (up to 0(a) corrections). 
Next, to determine Z^om in Eq. (13. 3p . we may use the conservation law 
Eq. (ETT51) . itself. 

The determinations of Z 2 Zg and c are somewhat correlated. We first 
note that, when the energy-momentum tensor T^ v (x) (13. 4p is inserted in a 



of 0^(x 



-'Ref. [65] is a pioneering work on the trace anomaly in lattice gauge theory. 
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physical amplitude, the last term of Eq. ( 13. 4p just gives rise to an additive 
constant (13. 7p that is identical for any physical amplitudes; the constant c 
thus corresponds to a choice of the origin of the energy. If we consider the 
difference in expectation values of T^^x) in two physical states, therefore, 
the contribution of the last term of Eq. (13 .4p cancels out and it is the same 
as the difference in expectation values of (1/2) [0^ u (x) + UfJi (x)]. Noting 
this fact, one may determine Z 2 Zs, the absolute normalization of 6> MJ ,(x) 
in Eq. (13. 3p . by setting the difference of expectation values of the "energy 
operator" — a 3 J2x @oo(%) i n two different physical states to a certain pre- 
scribed value. 

Although the value of c influences neither on the conservation law nor on 
the difference in expectation values of the energy-momentum tensor, there 
exists a physically natural choice of c in the present supersymmetric system. 
That is, we may require that the expectation value of the energy density to 
vanish 

(Too(x)) = 0, (3.21) 

when periodic boundary conditions are imposed on all the fields. This re- 
quirement fixes 

a 1 



2(N 2 — 1) ^ 00 ^ a '^P eriodic boundary conditions ' (3.22) 

because of Eq. (13. 7p . Eq. ( I3.2ip states that the derivative of the super- 



symmetric partition function (i.e., the Witten index [67|) with respect to 
the temporal size of the system vanishes; this holds only if a choice of the 
origin of the energy is consistent with the SUSY algebra. In other words, 
Eq. (I3.2ip is required for the spatial integral —a 3 J2x Too 0*0 t° be the energy 
operator appearing in the right-hand side of the SUSY algebra. It is thus 
a natural requirement from the perspective of SUSY0 This completes our 



n We emphasize that Eq. (|3.2ip should hold even in a theory in which SUSY is spon- 
taneously broken (recall that the Witten index is independent of the temporal size of the 
system even in such a case) . This same idea was adopted to find an expression of the en ergy 
density in a lattice formulation of the two- dimensional (2D) TV = (2, 2) SYM [HI M, MM- 
For this 2D system, fortunately, there exist lattice formulations that possess one exact 
fermionic symmetry [Til I72I l73l . [73 ] and, employing this symmetry, one can define the en- 
ergy density operator that shows the zero-point energy consistent with SUSY. This energy 
density operator on the lattice has been used 75[ to measure the vacuum energy density 
associated with a possible dynamical SUSY breaking in the above 2D system. 
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non-perturbative construction of a symmetric energy-momentum tensor. It 
is conserved in the quantum continuum limit and, by choosing c as Eq. (13. 22 p . 
the origin of the energy is consistent with SUSY. 

A remaining important issue is whether our symmetric energy-momentum 
tensor (13. 4p generates, through a renormalized WT relation, a correctly- 
normalized translation on renormalized fields. Although we do not go into 
this question in the present paper, assuming the existence of a renormalized 
translational invariant theory and considering the uniqueness of the con- 



served symmetric energy-momentum tensor in the continuum limit [48 



we 



believe that the answer is affirmative: Of course, further consideration is 
required on this point. 



4. Conclusion 



We have presented a non-perturbative construction of a symmetric energy- 
momentum tensor in the lattice formulation of 4D M = 1 SYM. Inspired by 
the relation (II .9p in the FZ supermultiplet, we defined the energy- momentum 
tensor by a renormalized modified super transformation of a renormalized 
SUSY current that generates a finite super transformation on renormalized 
operators (under the tuning of the gluino mass). Then, it can be shown that 
the lattice energy-momentum tensor is conserved in the quantum contin- 
uum limit. The resulting energy- momentum tensor may be used to measure 
physical quantities related with the energy-momentum tensor (such as the 
viscosity) by numerical simulations. 

A question naturally arises is then what will be resulted from the con- 
sideration on the another relation in the FZ multiplet, Eq. (II .8p . a relation 
between the U(1)a current and the SUSY current. It is clear that the clas- 
sical relation (I1.8P as it stands cannot hold in quantum theory, because the 
U(1)a current is not conserved owing to the axial anomaly, while we are 
believing that SUSY does not suffer from any anomaly. Ironically, it turns 
out that the situation for this seemingly simple relation (II. 8p is much more 
complicated than for (seemingly complicated) Eq. (jl.9P ; this is precisely be- 
cause of the axial anomaly. We hope to come back this problem in the near 
future. If the full structure of the FZ multiplet (with corrections by the 
anomaly) can be realized in a well-regularized framework such as the lattice, 
it will be a quite useful starting point to understand the so-called anomaly 



puzzle [76|, mi, lZ8|, [79|, |80|, |8l|, |82 



83, 84, 85, 86 
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Note added in proof 



Eq. ( I3.13P implies the following very suggestive relation: Noting that 
J2 y ev x duSviv) i n the first term of the right-hand side of Eq. (13.131) does not 
have the support at y = x, using the SUSY WT relation (I2.25j) again, we 
have 



« 4 E { d l°^)°) 



x IZ 



xev yen. 
1 d 



a A d£{x) 



A 5 + a£(x) 



Z 



1 d . . 



« 4 d^y) 



o 

(4.1) 



where T>q is a finite region that contains the operator O entirely. On the 
other hand, for the anti-symmetric part of 6> Mi ,(x), A^ v (x) in Eq. (13.161) . 



« 4 E ( d "A^)o) ^ 0, 

xex> 



(4.2) 



because Eq. (I3.17P shows that (A^ u (x)0) is proportional to the delta func- 
tion 5 4 (x — z) in the continuum limit, where z is any point in the support 
of O. Similarly, for the last term of Eq. (13 .4p . we have 



fl4 E (^<V tr + M)if>{x)] O) = 0. 



(4.3) 



Thus, combining the above three relations, 

« 4 E ( d %^)°) 
> -|(c -1 7.)^a 4 E ° 4 E 



a->0 



xeT> y<=T> x 



X ^ 



1 d 



Z 



1 g 



A 5 + a£(y) 



(4.4) 
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This may be regarded as the SUSY algebra in the present lattice framework 
if we assume that the lattice energy-momentum tensor T^ v generates the 
translation (on gauge invariant operators O) as 

« 4 E W%»( x )°) ^ - (W) . (4.5) 

Remark in proof 

A quite similar but somewhat different definition of a lattice energy- 
momentum tensor in 4D M = 1 SYM, which is superior in several aspects 



compared with the one in the present note, has been given in Ref. [87 . 
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